Dark halo densities, substructure, and the initial power spectrum 
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Although the currently favored cold dark matter plus cosmological constant model for structure formation 
assumes an n — 1 scale-invariant initial power spectrum, most inflation models produce at least mild deviations 
from n = 1. Because the lever arm from the CMB normalization to galaxy scales is long, even a small "tilt" can 
have important implications for galactic observations. Here we calculate the COBE-normalized power spectra 
for several well-motivated models of inflation and compute implications for the substructure content and central 
densities of galaxy halos. Using an analytic model, normalized against N-body simulations, we show that while 
halos in the standard (n — 1) model are overdense by a factor of ~ 6 compared to observations, several of our 
example inflation+LCDM models predict halo densities well within the range of observations, which prefer models 
with n ~ 0.85. We go on to use a semi-analytic model (also normalized against N-body simulations) to follow 
the merger histories of galaxy-sized halos and track the orbital decay, disruption, and evolution of the merging 
substructure. Models with n ~ 0.85 predict a factor of ~ 3 fewer subhalos at a fixed circular velocity than the 
standard n = 1 case. Although this level of reduction does not resolve the "dwarf satellite problem", it does 
imply that the level of feedback required to match the observed number of dwarfs is sensitive to the initial power 
spectrum. Finally, the fraction of galaxy-halo mass that is bound up in substructure is consistent with limits 
imposed by multiply imaged quasars for all models considered: f sa t > 0.01 even for an effective tilt of n ~ 0.8. We 
conclude that, at their current level, lensing constraints of this kind do not provide strong limits on the primordial 
power spectrum. 



1. Introduction 

The cold dark matter plus cosmological con- 
stant (LCDM) model of structure formation is 
the most successful and popular current theory 
for the origin of universal structure. The the- 
ory is highly predictive, and given that obser- 
vations constrain the universe to be nearly flat 
(with fi m = 1 — S1a ~ 0.3, and h « 0.7) a 
crucial unknown is the primordial spectrum of 
density fluctuations. It is usually assumed that 
the initial fluctuation spectrum is scale- invariant, 
with P(k) oc k n , n = 1. This choice is moti- 
vated by the fact that slow-roll inflation mod- 
els predict nearly scale-invariant spectra. How- 
ever, generic models of inflation do not predict 
primordial power spectra that are exactly scale- 
invariant and almost all models have some small, 
often non-negligible, "tilt" (n ^ 1) and "run- 
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ning" of the spectral index (dn/dhik ^ 0). In 
this proceeding we report on our work to model 
the small-scale implications of reasonable initial 
power spectra, and demonstrate that the sub- 
structure content and central densities of dark 
halos should be very sensitive to small devia- 
tions from scale-invariance. Our models are semi- 
analytic, but they are normalized and tested 
against N-body simulations. As we discuss be- 
low, our results may be important for interpret- 
ing galaxy rotation curves, dwarf galaxy counts, 
and substructure mass fractions as measured from 
multiply imaged quasars. 

In Zentner & Bullock ||| (hereafter ZB1 and 
ZB2) we chose several well-motivated and rep- 
resentative inflation models and calculated the 
implied power spectra to second order in slow- 
roll using the method of Stewart and Gong fll4|] . 
Adopting fi M = 0.3, fi L = 0.7, tt h h 2 = 0.020, and 
h = 0.72 (LCDM), we calculated P(k) using the 
results of Eisenstein & Hu |5f| and normalized to 
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COBE on the pivot scale, k* ~ 0.0023 h^Mpc. 
In the interest of brevity, we report on only three 
models here: 

• IPL: n = 0.94, dn/dlnk = -0.001, 
a 8 = 0.83 

• RM: n = 0.84, dn/d\nk = -0.004, 
a & = 0.65 

• n = 1: n = 1.00, dn/d\nk = 0.000, 
as = 0.95. 

The tilt and running are evaluated at k — /c* and 
as is the rms fluctuation amplitude within spheres 
of 8/i~ 1 Mpc implied by the COBE normalization. 
Observational constraints on as have historically 
favored a value near unity; however, some recent 
estimates indicate surprisingly low values of as 
||. Roughly speaking, 0.60 < a s < 1.2 spans the 
range of recently reported values (for il m ~ 0.3). 

2. Central Densities 

The theory of CDM has survived remarkably 
well since its creation roughly twenty years ago 
but lingering for nearly half this time has 
been a possible problem with galaxy central den- 
sities H , which seem to be much less dense than 
the standard CDM model typically predicts. One 
part of this problem — the cusp vs. core question 
— concerns the precise inner slope of galaxy rota- 
tion curves and whether they match the "cuspy" 
predictions of CDM theory. This particular ques- 
tion is highly debated, and some argue that ob- 
servational errors and/or theoretical uncertain- 
ties make the comparisons inconclusive. But even 
among those who believe there may be no prob- 
lem with the shapes of rotation curves (the au- 
thors included), it is almost universally agreed 
upon that the data do prefer halos that are less 
centrally concentrated than typical halos in the 
standard (n = 1) LCDM model. 

One way to see the problem is by comparing 
rotation curve data to theory using the following, 
non-parametric measure of the density m : 
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Figure 1. The central density parameter (see 
text) as a function of maximum circular veloc- 
ity for several dark matter-dominated galaxies 
(points). The solid line corresponds to the pre- 
dicted densities for standard n = 1 LCDM halos, 
and the short and long-dashed lines correspond 
to our expectations for two different inflation- 
derived input power spectra. 



This is the mean dark matter density 
within the radius ?V/2 where the galaxy 
rotation curve falls to half its maximum, 
Vmax- In practical units, Ay/2 — 5 x 
10 5 (Knax/100 kms-^trv/a/rr 1 kpc)" 2 . 

In Figure 1, we plot A y/ 2 measured from 
several dark matter-dominated galaxy rotation 
curves (data from fllql , see ZB1 for details) and 
compare to predictions for typical halo densities 
in our models with different power spectra (lines) . 
Notice that the n — 1 model overpredicts the 
galactic densities by a factor of ~ 6, while our 
inflation-motivated models with n < 1 do much 
better. 

The predictions were made using the semi- 
analytic model of Bullock et al. which was 
shown to reproduce the halo densities measured 
in N-body simulations for n = 1 and 0.9 LCDM 
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power spectra as well as a variety of power-law 
models (see ZB1 for an extensive discussion). The 
model is similar to that originally proposed in 
11 1 (see also @), although modified in order 
to reproduce the redshift dependence seen in N- 
body halos. Qualitatively, it embodies the fact 
that halo densities reflect the density of the uni- 
verse at their characteristic collapse redshift. Ha- 
los in models with less small-scale power collapse 
later, when the universe is less dense, and there- 
fore have more diffuse cores. Interestingly, even 
the "IPL" case, with only a moderate tilt and 
mild running, does significantly better than the 
n = 1 standard model. If we adopt the crite- 
rion that a "good" fit to the data has Ay/ 2 
3 x 10 5 at y max 



1 then the data prefer 
0.85, assuming that 



100 kms" 
n(k*) + 6.96dn(k*) /dhik 
the primordial spectrum is the only element of 
the solution. 

3. Substructure 

As mentioned briefly in the previous section, 
the input power spectrum greatly affects the typ- 
ical merger history of galaxy-sized halos. Ha- 
los formed in models with less power on galactic 
scales typically form later, and accrete a larger 
fraction of their substructure at later epochs. 
Similarly, the "subhalos" that are accreted also 
form later, are less dense, and therefore are more 
easily tidally disrupted. In order to explore how 
different initial power spectra affect the substruc- 
ture content of halos quantitatively, we developed 
a semi-analytic model that tracks halo accretion 
histories using the extended Press-Schechter for- 
malism [|l3| and follows subhalo orbital evolution 
in a semi-analytic fashion. 

Described in detail in ZB2, our model ex- 
pands on the calculations presented in Bullock, 
Kravtsov & Weinberg [|| and includes some mod- 
ifications similar to those advocated by Taylor 
and Babul |16| . Briefly, after each merger event, 
we follow the subhalo orbital evolution as it de- 
cays due to dynamical friction and loses mass via 
tidal interactions with the host halo. Dynami- 
cal friction dominates for large subhalos, which 
typically merge with the central object soon after 
accretion. Low mass substructures are efficiently 
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Figure 2. Number of subhalos as a function of cir- 
cular velocity. The functions are derived for stan- 
dard (n = 1) LCDM host halos of mass Afyir — 
10 12 ft.~ 1 M©. The long-dashed line shows all ac- 
creted subhalos. The short-dashed line shows 
subhalos that were tidally destroyed, and the dot- 
ted line those that merged with the central ob- 
ject. Subhalos that exist at z — are shown by 
the thick solid line. Error bars correspond to the 
rms scatter. Subhalo counts from two different 
N-body simulations are shown by the thin solid 
lines. Both agree reasonably well with our model 
expectations. 



stripped by the host potential, and a large frac- 
tion of them become completely tidally destroyed. 
The results of our model for an ensemble of Milky 
Way-like host halos in an n — 1 LCDM cosmology 
are displayed in Figure 2. Shown as a function 
of V„iax are the total number of objects that have 
been accreted compared to the number that even- 
tually merge with the central object, the number 
that are tidally destroyed, and the number that 
survive as distinct substructures at z = 0. It is 
clear that our model for the number of surviving 
structures is in good agreement with the results of 
N-body simulations (thin solid lines) B. We use 
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Figure 3. The final subhalo abundances for our 
three models compared to two estimates of the 
local population of dwarf galaxies (points). 



this model to estimate the sensitivity of the sub- 
structure population at z = to the primordial 
power spectrum. 

Shown in Figure 3 are the cumulative subhalo 
velocity functions for each of our models com- 
pared to two different estimates of dwarf galaxy 
counts within ~ 200kpc of the Milky Way and 
M31 |J|. That n = 1 LCDM overpredicts the 
number of small satellites by roughly an order of 
magnitude is expected. This is the well-known 
dwarf satellite problem (e.g. ||). Interestingly, 
the number of subhalos is reduced by as much as a 
factor of ~ 3 for the inflation-derived power spec- 
tra that we consider. Although all of these cases 
would require some degree of feedback to match 
dwarf counts (e.g., ionization ||), the amount of 
feedback needed will depend sensitively on the in- 
put power spectrum. 

Related to the dwarf satellite problem are re- 
cent efforts to measure the substructure content 
in halos using the flux ratios in multiply imaged 
quasars po[ . For example, Dalai and Kochanek 
used this technique to constrain the fraction of 



Figure 4. The fraction of mass in a galaxy-host 
halo that is contained in substructure less massive 
than M. The box roughly represents the lensing 
constraints of Dalai and Kochanck||To[|. 



the host halo mass bound up in substructure to 
be 0.006 < / sat < 0.07 (90 % confidence). They 
went on to use this result to try to limit the 
tilt of the primordial spectrum using a simpli- 
fied model for the abundance of substructure and 
obtained n > 0.94 at 95% confidence. Our re- 
sults on substructure differ from those of Dalai 
and Kochanek, as illustrated in Figure 4. We find 
that for a host halo of the relevant mass, the total 
mass fraction in subhalos is typically larger than 
the their lower limit (/ sa t > 0.006) even for our 
most tilted primordial spectrum, n ~ 0.84. Since 
this model is on the edge of what is acceptable by 
cr 8 constraints, it seems that strong lensing does 
not yet serve as a significant probe of the pri- 
mordial power spectrum. (Using a version of our 
model, Dalai and Kochanek are in the process of 
revising their limits, and now find results more in 
line with those reported here H.) 
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